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Explicit expressions for restricted partition function W(s,d m ) and its quasiperiodic compo- 
nents Wj(s, d m ) (called Sylvester waves) for a set of positive integers d m = {d\, da, ■ ■ ■ , d m } are 
derived. The formulas are represented in a form of a finite sum over Bernoulli polynomials of 
' higher order with periodic coefficients. 

1 Introduction 
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Abstract 



The problem of partitions of positive integers has long history started from the work of Euler who 
laid a foundation of the theory of partitions p^, introducing the idea of generating functions. Many 
prominent mathematicians contributed to the development of the theory using this approach. 

J.J. Sylvester provided a new insight and made a remarkable progress in this field. He found 
jl()l 111] the procedure for computation of a restricted partition functions, and described symmetry 
properties of such functions. The restricted partition function W(s,d m ) = W(s, {d\, di, • . . ,d m }) 
is a number of partitions of s into positive integers {d\,d2, ■ ■ ■ , d m }, each not greater than s. The 
generating function for W(s,d m ) has a form 

a" 
> 

X 

' where W(s, d m ) satisfies the basic recursive relation 



m 1 oo 

i? ^ dm )=n 1 ^-=E^ dm )^' « 



W(s, d m ) - W(s - d m ,d m ) = W(s, d m - L ) . (2) 

Sylvester also proved the statement about splitting of the partition function into periodic and non- 
periodic parts and showed that the restricted partition function may be presented as a sum of 
"waves", which we call the Sylvester waves 

W(s,d m ) = ^2w j (s,d m ), (3) 

where summation runs over all distinct factors in the set d m . The wave Wj(s,d m ) is a quasipoly- 
nomial in s closely related to prime roots pj of unit. Namely, Sylvester showed in that the wave 
Wj(s,d m ) is a coefficient of t^ 1 in the series expansion in ascending powers of t of 

p~ s e st 

F J st ) = Y p r . (4) 
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The summation is made over all prime roots of unit pj = exp(2irin/j) for n relatively prime to j 
(including unity) and smaller than j. This result is just a recipe for calculation of the partition 
function and it does not provide explicit formula. 

Using the Sylvester recipe we found in [2] an explicit formula for the Sylvester wave Wj(s,d m ) 
in a form of finite sum of the Bernoulli polynomials of higher order [21 [7j multiplied by a periodic 
function of integer period j. The periodic factor is expressed through the Eulerian polynomials of 
higher order [I]. In this note we show that it is possible to express the partition function through 
the Bernoulli polynomials only. 

A special symbolic technique is developed in the theory of polynomials of higher order, which 
significantly simplifies computations performed with these polynomials. A short description of this 
technique required for better understanding of this paper is given in the Appendix. 



2 Polynomial part of partition function and 
Bernoulli polynomials 

Consider a polynomial part of the partition function corresponding to the wave VKiXs, d m ). It may 
be found as a residue of the generator 

Fl(s,t) = T-rm n Zrf^V ■ ( 5 ) 

n^iC 1 — e i ) 

Recalling the generating function for the Bernoulli polynomials of higher order [2]: 

st+mytm , 00 n 

e 1 lli=l a i _ ST^ p>(m)(„\ A m\ l 

nr=i(^-i) to n 

and a transformation rule 

B^\s\-d m )=B^\s + Y,dM 

i=i 

we obtain the relation 

„st 1 00 4.n—m 

— ^2B^(s + s m \d m )—, (6) 
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where 

m 

i=l i=l 

It is immediately seen from (JHJ) that the coefficient of 1/t in © is given by the term with 
n = m — 1 

W ^ dm ) = (m _\ Vv B ™-i( s + s ™l d ' m ) • ( 7 ) 

yiil 1). 7T m 

The polynomial part also admits a symbolic form frequently used in theory of higher order polyno- 
mials 

/ m \ m-1 



i=l 
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where after expansion powers rj of l B are converted into orders of the Bernoulli numbers % B Ti =>■ B ri . 
It is easy to recognize in the explicit formula reported recently in which was obtained by a 
straightforward computation of the complex residue of the generator (J5J) . 
Note that basic recursive relation for the Bernoulli polynomials (Zj 

4 m) (s + d m \d m ) - B^\s\d m ) = n ( i mJ B^- 1) ( S |d m - 1 ) 

naturally leads to the basic recursive relation for the polynomial part of the partition function: 

W l (s,d m )-W 1 (s-d m ,d m ) = W 1 (s,d m - 1 ) , 

which coincides with Q. This indicates that the Bernoulli polynomials of higher order represent a 
natural basis for expansion of the partition function and its waves. 

3 Sylvester waves and Eulerian polynomials 

Frobenius [Hj studied in great detail the so-called Eulerian polynomials H n (s,p) satisfying the gen- 
erating function 

H n=0 

which reduces to definition of the Euler polynomials at fixed value of the parameter p 

E n (s) = H n (s,-l). 

The polynomials H n (p) = H n (0,p) satisfy the symbolic recursion (Hq(p) = 1) 

pH n (p) = (H(p) + 1)", n > 0. (8) 
The generalization to higher orders is straightforward 

p sty\ m (-\ _ n di\ °° f n 

Y\i=i(e r -P l ) ^ n\ 

where the corresponding recursive relation for H^ n \s, p\d m ) has the form 

Ht\s + d m ,p\d m ) - p d -H^(s,p\d m ) = (l - HlT-V&pld™- 1 ) . 

The Eulerian polynomials of higher order Hn \s, p\d m ) introduced by L. Carlitz in [I] can be 
defined through the symbolic formula 

/ m \ n 
Ht\s,p\d m )= (s + J^di'Hip^)) , 



i=l 



where H n (p) computed from the relation 

r n=0 

or using the recursion (JSJ). Using the polynomials Hn \s, p\d m ) we can compute the Sylvester wave 
of arbitrary period. 
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In order to compute the Sylvester wave with period j > 1 we note that the summand in the 
expression (JIJ) can be rewritten as a product 



,,st 



n£i(l-e-^) im* J+1 (l-pfe-**)' 



(9) 



where the elements in d m are sorted in a way that j is a divisor for first kj elements (we say that j 
has weight kj), and the rest elements in the set are not divisible by j. 

Consider the j-periodic Sylvester wave Wj(s,d m ), and rewrite the summand in Q as double 
infinite sum 



— j- ^ D n {S + S k] |d 

lli=fc,-+l(l ~ Pj) n=0 



j2Hl m -^( Sm -s k ., Pj \d^y- 



t-j ni=ij + n- i'j y n=U 1=0 

The coefficient of 1/t in the above series is found for n + I = kj — 1, so that we obtain a finite sum: 

-s 
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This expression may be rewritten as a symbolic power : 

— S I 771 

W^d") = , k ._\ y E U m odt) ls + S m + J2 d i iB + J2 d * ^ 



(11) 



which is equal to 



Wj(s,d r 
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(12) 



where 



^V|d m ]=ffM(0,p|d' 



E^^(p d 



i=l 



are the Eulerian numbers of higher order and it is assumed that 

#< O) |>|0] = 1, /4°)[p|0]=O, n>0. 

It should be underlined that the presentation of the Sylvester wave as a finite sum of the Bernoulli 
polynomials with periodic coefficients ()12|) is not unique. The symbolic formula (|llj) can be cast 
into a sum of the Eulerian polynomials 



kn-l 



Wj{s,d r ' 



(kj - 1)! 7Tfc 
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J n=0 
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(13) 
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where 

#M[ d m] =jB M( (d m ) 

are the Bernoulli numbers of higher order. 

It should be noted that the formulae (|12|) and (jl.3[) require summation over all prime roots of unit, 
and though it is simpler than the Sylvester recipe using @, it cannot be considered a completely 
explicit formula. 



4 Reduction of Sylvester waves to Bernoulli polynomials 

A relation between the Eulerian and Bernoulli numbers and polynomials of higher order was estab- 
lished in [1], it may be written as follows: 

(m-l-n)! n m Pj m 1 rr{m-k 3 ) 



(kj-l-ny.n^pp Il? =kj+I (l-Pp "' 



(m j—1 \ m 

i=kj+lri=0 J i=kj+l 

Using this relation we convert the inner sum in (|13|) into 

y — — — ^H { ^i{s + Sm , P ,id m - fc o 

Pj Hi=fc :( +iv i Pj > 



j 
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(kj - 1 - n)\ 7T fe , 
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(m - 1 - n)\ ir m 

Pi 

m j—1 \ m 

=0 J 

(kj - 1 - n)\ ix kj 



-(m—kj) 



(m — 1 — n)\ Tr n 



m J—1 
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Here ^ j{s) denotes a j-periodic prime radical circulator introduced in [5] 



*j(s)=J2p. 

For prime j it is given by 



3' 

Pi 



<•«> 

where </>(j) and denote Euler totient and Mobius functions. Considering j as a product of 
powers of distinct prime factors 

3 =n^ fe ' 
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one may easily check that for integer values of s 



%(*)=npr%j^j> (17) 



where VPfc(s) = for non-integer values of s. 

/ xi >i I tl vl 1 III III ;l Hi] K 1 \ 



It is convenient to introduce a j -modified set of summands d^ 1 defined as union of subset d fc ^ of 



summands divisible by j and the remaining part d m k i multiplied by j 

so that d™ is divisible by j. Substitution of (|15|) into ()13j) with extension of the outer summation 
up to m produces the formula for computation of the j-periodic Sylvester wave: 

m j-1 



m 

x B^} 1 (s + s m +J2 ^r i |dJ l )* j ( S + E diin + l)). (18) 

i=kj+l i=kj+l 

The derivation of (|18|) implies that all terms containing s in powers larger than kj — 1 are identically 
equal to zero. The polynomial part of the partition function W\{s , d™) for the j-modified set of 
summands reads: 

Wl(s,d7) = {m _ lV \ ^J ^l (« + *m + - 1) E *|dj»). 

This formula gives rise to the representation of the j-periodic Sylvester wave Wj(s, d m ) through the 
linear combination of the polynomial part of the j-modified set of summand d™ multiplied by the 
j-periodic functions ^fj-. 

(m j \ m m 

II E )Wi(s+ E d i (r i -j),d?)y j (s+ E diin-j)), 
i=kj+l Ti=\ J i=kj+l i=kj+l 

which is written also as 

(m j— 1 \ m m 

H e wi( S -E^ d f)%( s -E^)- ( 19 ) 
i=kj+lri=Q/ i=kj+l i=kj+l 

The last formula shows that each Sylvester wave is expressed as a linear superposition of the polyno- 
mial parts of the modified set of summands multiplied by the corresponding prime circulator. Thus, 
the formulae (|16ll8j) with the Sylvester splitting formula © provide a complete explicit solution of 
the restricted partitions problem. 



Appendix 

The symbolic technique for manipulating sums with binomial coefficients by expanding polynomials 
and then replacing powers by subscripts was developed in nineteenth century by Blissard. It has 
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been known as symbolic notation and the classical umbral calculus 8 . An example of this notation 
is also found in |2j in section devoted to the Bernoulli polynomials B^{x). 
The well-known formulas 

B n (x + y) = £ ( n ^B k (x)y n -\ B n (x) = £ (f\ B k^ n ~\ 

k=0 ^ ' k=0 ^ ' 

are written symbolically as 

B n (x + y) = (B(x) + y) n , B n (x) = (B + x) n . 

After the expansion the exponents of B(x) and B are converted into the orders of the Bernoulli 
polynomial and the Bernoulli number, respectively: 

[B(x)] k B k (x), B k ^B k . 

We use this notation in its extended version suggested in J7| in order to make derivation more clear 
and intelligible. Norlund introduced the Bernoulli polynomials of higher order defined through the 
recursion 

Bt\x\d™) = £ ( n \d k B k {Q)Btl l \x\d m - 1 \ 
k=o ^ ' 

starting from Bn {x\d\) = d™B n (-^). In symbolic notation it takes form 

B^(i)=(4B(0)+flM(i)) tt I 
and recursively reduces to more symmetric expression 

(m \ n 

where each [ l B(0)] k is converted into B k (0). 
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